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VIBRATION  ANALYSIS  OF 
CLUSTERED  LAUNCH  VEHICLES 


SUMMARY 


Irhe  Saturn  I  vehicle  represents  a  complex  ensemble  of  nine  beams  joined  together. 
A  first  vibration  analysis  based  on  an  equivalent  single  beam  did  not  agree  satisfactorily 
with  the  test  results  of  the  full  scale  test  vehicle.  Comparison  of  analysis  and  test  re¬ 
sults  rather  dictates  that  the  coupling  effects  of  the  nine  single  beam  vibrations  can  not 
be  neglected.  Hence,  the  analysis  developed  in  this  paper  is  based  on  a  model  consisting 
of  nine  clustered  beams J  However,  the  complexity  of  the  system  calls  for  many  simpli¬ 
fying  assumptions. 

Irhe  basic  idea  of  the  analysis  consists  in  coupling  of  the  single  beam  vibrations 
such  that  the  connections  of  the  beams  will  be  preserved.  The  method  applied  to  the  un¬ 
coupled  single  beams  is  the  well  known  transfer  matrix  method.  The  mentioned  coupling 
of  the  single  beam  vibrations  leads  to  a  linear  homogeneous  system  of  compatibility 
equations  which  must  be  solved  during  each  transfer.  Realization  of  the  boundary  con¬ 
ditions  leads  to  a  frequency  curve  indicating  the  coupled  frequencies  of  the  system. 

Finally  the  coupled  mode  shapes  can  be  determin^J 


I^The  analysis  is  programmed  on  IBM-709-4.  Because  it  consists  of  simple 
matrix  operations  the  computer  time  is  short,  j 


The  results  of  an  analysis  of  the  Saturn  I  vehicle  compare  well  with  test  results,  | 


SECTION  1.  INTRODUCTION 


The  Saturn  vehicle  as  shown  by  Figure  1  consists  of  payload,  S-IV  stage  and  S-I 
stage.  In  the  S-I  stage,  eight  clustered  tanks  (four  fuel  and  four  lox)  are  attached  at 
each  end  to  the  center  barrel  (lox  tank) .  Then,  the  Saturn  vehicle  represents  a  highly 
complex  ensemble  of  nine  beams  joined  together. 

The  results  of  a  first  vibration  analysis  based  on  an  equivalent  single  beam  did 
not  agree  satisfactorily  with  the  test  results  of  the  full  scale  test  vehicle  (Ref.  1).  The 
comparison  of  analysis  and  test  results  dictates  that  coupling  can  not  be  neglected;  indeed 
the  vibrations  of  the  nine  single  beams  must  be  coupled. 


Subsequently,  several  proposals  concerning  the  vibration  analysis  of  clustered 
booster  vehicles  have  been  made.  Some  of  them  are  suggested  by  References  2,  3,  and  4. 
The  analysis  shown  herein  is  that  done  by  the  Dynamics  and  Loads  Branch  (Structures 
Division  of  the  Propulsion  and  Vehicle  Engineering  Laboratory)  George  C.  Marshall 
Space  Flight  Center. 

This  analysis  is  of  a  model  outlined  in  Section  III  of  this  report.  Complexity  of 
the  Saturn  vehicle  calls  for  many  simplifying  assumptions.  But,  remember  the  funda¬ 
mental  premise:  the  basic  configuration  of  the  model  is  a  central  beam  with  eight 
smaller  beams  clustered  about  it.  From  this  premise,  simplifications  can  be  made 
concerning  the  mode  of  vibration,  and  the  support  conditions  of  the  outer  beams  and 
single  beam  models.  Assumptions  about  the  mode  of  vibration  are: 

a.  The  center  beam  is  restricted  to  plane  vibrations  in  a  symmetry  plane  of  the 
vehicle  through  its  long  axis. 

b.  The  outer  beam  vibrations  occur  symmetrically  to  this  plane. 

Nevertheless,  each  outer  beam  may  have  two  planes  of  vibrations.  The  outer  beam 
supports  are  assumed  equivalent  to  either  simple  supports,  or  to  simple  supports  having 
additional  bending  restraints.  The  single  beams,  finally,  are  modeled  by  lumped  mass 
systems  based  on  Timoshenko^  s  beam  theory. 

The  essential  part  of  the  analysis  (Section  IV)  deals  with  the  coupling  of  the 
single  beam  vibrations.  The  beams  do  not  vibrate  independently — they  are  connected. 
Coupling  means  adjustment  of  outer  beam  and  center  beam  vibrations  such  that  the 
connections  of  the  beams  will  be  preserved.  In  other  words,  certain  conditions  of  com¬ 
patibility  must  be  satisfied.  These  circumstances  can  be  expressed  mathematically  by 
a  linear  inhomogeneous  system  of  equations  to  be  solved. 

The. fundamental  method  applied  in  this  analysis  is  the  transfer  matrix  method. 

This  method  is  well  known  (Ref.  5)  but,  for  completeness  a  short  explanation  is  given  in 
the  appendix.  The  transfer  matrix  method  is  based  on  the  assumption  that  the  beam 
masses  and  mass  moments  of  inertia  are  lumped  at  stations  of  the  beam  axes.  The 
state  at  each  station  is  determined  by  a  state  vector  which  contains  deflection,  rotation 
angle,  moment  and  shear  at  this  station  as  components.  The  state  vectors  of  successive 
stations  are  linked  by  transfer  matrices.  Such  a  transfer  matrix  is  the  product  of  an 
inertia  and  a  stiffness  matrix.  In  this  way  the  transfer  from  one  station  to  any  other 
station  of  the  beam  is  determined  by  a  matrix  which  is  composed  of  inertia  and  stiffness 
matrices.  This  applies  especially  to  the  end  stations  of  the  beam.  These  beam  transfer 
matrices  of  the  center  and  the  outer  beams  are  defined  in  Table  2,  From  the  elements  of 
these  matrices  the  above  mentioned  compatibility  equations  can  be  formed  and  finally 
solved.  Then  the  transfer  of  the  complete  system  can  be  performed  (Table  3).  Now, 
realization  of  the  boundary  conditions  leads  to  the  eigenvalue  equation.  Evaluation  of  this 
equation  and  determination  of  the  mode  shapes  form  the  final  steps  of  the  analysis. 


2 


The  whole  procedure  outlined  above  will  be  done  by  application  of  simple  matrix 
algebra.  Formation,  multiplication,  inversion  of  low  order  matrices  are  the  only  opera¬ 
tions  performed.  The  matrices  used  and  their  dimensions  are  compiled  in  Table  4.  The 
analysis  is  programmed  on  IBM-709-4.  Because  it  consists  of  the  mentioned  simple 
matrix  operations,  the  computer  time  is  very  short.  To  obtain  one  mode  shape  and  fre¬ 
quency,  approximately  three  minutes  are  necessary. 

The  results  of  an  analysis  of  the  Saturn  I,  SA-5  vehicle  are  presented  and  compare 
well  with  the  test  results.  The  mode  shape  characteristics  of  these  complex  shapes  in¬ 
dicate  the  proper  coupling  model  is  used  in  the  analysis. 


SECTION  II.  MODEL 

Before  beginning  discussion  of  the  model  used,  a  short  description  of  the  vehicle 
will  be  given  (Fig.  1) .  The  vehicle  consists  of  payload,  second  stage  and  first  stage.  In 
the  first  stage,  eight  tanks  are  clustered  symmetrically  about  a  center  tank.  To  connect 
outer  tanks  with  the  center  tank,  a  spider  beam  at  the  upper  end  and  the  outrigger  assembly 
on  the  lower  end  are  used  (Fig.  2) .  In  addition,  spider  beam  and  outrigger  assembly 
provide  interconnection  between  stages  and  supports  of  the  engines,  respectively.  The 
eight  outer  tanks  consist  of  four  lox  tanks  and  four  fuel  tanks .  The  center  tank  is  a 
lox  tank. 

The  basic  configuration  of  the  model  is,  therefore,  a  center  beam  with  eight 
smaller  beams  clustered  symmetrically  about  it.  Spider  and  outrigger  structure  provide 
supports  of  the  outer  beams. 

The  preparation  of  a  vibration  analysis  of  this  system  represents  a  difficult  pro¬ 
blem.  The  elastic,  vibrating  structure  consists  of  a  three-dimensional  ensemble  of  non- 
uniform  beams  connected  by  the  complex  spider  and  outrigger  structure.  Obviously,  the 
analysis  must  be  based  on  a  simplified  model.  On  the  other  hand,  analytic  results  based 
on  too-simple  models  do  not  agree  satisfactorily  with  the  test  results  of  the  full  scale 
test  vehicle.  Studies  made,  using  different  models,  indicate  that  realistic  assumptions 
concerning  the  outer  beam  supports  are  significant.  Thus,  to  build  up  a  mathematical 
model,  two  basic  principles  should  be  considered; 

a.  For  analytic  reasons  a  simplified  model  must  be  used. 

b.  To  represent  satisfactorily  the  dynamic  behavior  of  the  actual  vehicle,  the 
outer  beam  supports  of  the  spider  and  tail  sections  must  be  modeled  as  closely  as  possible. 

Some  of  the  simplifying  assumptions  are  usual  and  well  known: 

a.  The  nine  beams  of  the  vehicle  model  will  be  considered  as  Timoshenko  beams. 
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b.  The  beam  masses  will  be  lumped  at  stations  of  the  beam  axes  such  that  the 
stiffness  between  stations  may  be  assumed  as  constant. 

c.  The  liquid  within  the  tanks  will  be  assumed  as  solid — having  mass,  but  no 
mass  moment  of  inertia.  The  masses  will  be  lumped  like  the  other  solid  masses  on  the 
beam  axes. 

Another  assumption  deals  with  the  vibration  itself.  Obviously  a  three-dimensional 
cluster,  as  in  the  Saturn  vehicle,  may  vibrate  in  different  modes  depending  on  the  nature 
of  the  excitation.  For  simplification,  the  analysis  proposed  is  restricted  to  vibrations 
excited  in  a  plane  of  symmetry.  Cases  (a)  and  (b)  of  Figure  3  show  the  positions  of 
such  planes  within  the  cluster.  V  Then,  th^  center  beam  vibration  occurs  in  the  symmetry 
plane.  The  outer  beam  vibrates  symmetrically  to  this  plane.  Although  the  center  beam 
vibration  is  in  one  plane,  the  outer  beams  must  be  expected  to  vibrate  in  two  planes. 

The  analysis  under  discussion  is  based  on  case  (a)  of  Figure  3.  Yet,  the  analysis 
can  be  extended  to  case  (b) .  On  the. other  hand,  test  results  of  the  full  scale  test  vehicle 
do  not  indicate  considerable  deviations  of  mode  shapes  and  frequencies  of  the  cases  (a), 
(b) ,  and  (c)  shown  by  Figure  3. 

The  two  planes  of  vibration  of  the  outer  beams  can  be  assumed  as  radial  and 
tangential  planes  (Fig.  4a).  To  distinguish  among  the  different  vibration  planes,  these 
planes  will  be  numbered  (Fig.  4a).  For  reasons  of  symmetry,  restriction  to  one  side 
of  the  vehicle  is  possible.  The  planes  of  radial  vibrations  are: 

(1)  (2)  (4)  (6)  (8) 

while  tangential  vibrations  occur  in  the  planes: 

(3)  (5)  (7) 

To  measure  deflections  and  rotation  angles,  coordinate  systems  must  be  assigned 
to  each  vibration  plane  (Fig.  4b) .  Hence  the  outer  beam  vibration  planes  may  be 
characterized  by  the  angles  a,  =  1, 2,  ...  (8)  between  the  y-axes  of  the  center  beam  and 
outer  beam  vibration  planes.  Since  only  the  cosine  of  these  angles  is  important,  the 
angles  may  be  measured  in  arbitrary  directions  (Fig.  4a,  Table  1). 


1.  It  must  be  mentioned  that  these  planes  do  not  represent  perfect  symmetry  planes 
since  the  arrangement  of  the  outriggers  is  not  symmetric  to  these  planes  (Figs.  4  and  5) . 
However,  the  differences  between  both  types  of  outriggers  are  so  small — especially  with 
regard  to  stiffness — that  the  outriggers  can  be  considered  as  equal. 
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Now,  the  most  important  assumptions  have  to  do  with  the  outer  beam  supports 

located  on  spider  arms  and  outriggers  (Fig.  2).  Spider  arms  and  outriggers  will  be 

regarded  as  massless  beams  rigidly  connected  with  the  center  beam  at  two  branching 

points:  Branching  point  1  in  the  spider  section,  branching  point  2  in  the  tail  section. 

The  flexibility  of  spider  arms,  and  outriggers  can  be  determined  by  the  influence  co*- 

efficients  6  ,  6  respectively  of  the  support  points.  The  6's  are  the  deflections  per 

SP  0 

unit  load  in  direction  of  the  tank  axes.  In  view  of  the  outer  beam  connections  shown  by 
Figure  2,  the  following  assumptions  can  be  made. 


Spider  section 


In  the  radial  planes  of  vibration  (1),  (2),  (4),  (6),  (8),  and  in  the  tangential  vibration 
plane  (5)  the  outer  beams  behave  like  simple  supported  beams. 


Hence,  the  boundary  conditions  are 


(1) 


where  are  deflection  and  moment  of  the  i^^  vibration  at  the  spider  supports 

and  yi  is  the  center  beam  deflection  of  branching  point  1.  Hence  y^  costt^  represents  the 

support  deflection  in  vibration  plane  (i)  of  the  spider  section  caused  by  the  center  beam 
motion  of  branching  point  1. 

In  the  tangential  planes  of  vibration  (3)  and  (7)  the  outer  beams  behave  like 
supported  beams  having  additional  bending  restraints  ( spring  constant  x )  • 

This  leads  to  the  following  boundary  conditions: 


(2) 


where  <P-^^\  are  the  rotation  angles  of  the  i^^  outer  beam  and  the  center  beam  at 
branching  point  1  respectively. 
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Figure  5  shows  a  schematic  sketch  of  the  support  conditions  of  the  tangential 
vibration  planes.  For  convenience,  the  outer  beam  axes  (in  reality,  situated  on  a  cylinder 
around  the  center  beam)  are  rolled  up  in  the  drawing  plane.  One  concludes  that  the 
spring  stiffness 

d  2 
o 


where  d^  represents  the  outer  tank  diameter. 

Since  bending  of  the  center  beam  causes  rotation  of  the  spider  beam  plane  about 

(3)  (7) 

the  angle  the  angles  cp^  and  cp^  are  not  the  effective  angles  acting  on  the 
springs.  Rather  the  effective  angles  are  given  by  the  parenthetical  expression  of  the 
last  two  equations  of  the  system  2. 

Tail  Section 


In  the  radial  planes  of  vibration  (1),  (2),  (4),  (6),  (8)jthe  outer  beams  can  be  con¬ 
sidered  again  as  simple  supported  beams.  Hence,  the  boundary  conditions  of  these 
planes  may  be  written: 


(i) 

y2  =  y2  cos  a 

M,<‘>  =  0  1  =  1,  2,  4,  6,  8 


(3) 


where  y2^^^  >  represent  deflection  and  moment  of  the  i^^  outer  beam  vibration  at 

the  tail  support  and  y2  is  the  center  beam  deflection  of  branching  point  2. 


The  beams  vibrating  in  the  tangential  planes  (3),  (5),  (7), may  be  regarded  as 
supported  beams  having  additional  bending  restraints.  So,  with  regard  to  the  deflections 
it  follows 


y2 


(i) 


y2  cos  ;  i  =  3,  5,  7 


(4) 


To  obtain  the  remaining  boundary  conditions  of  the  vibration  planes  (3) ,  (5) ,  (7) , 
Figure  5  may  be  considered.  It  follows  easily  that 
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(5) 


<P2 

^2 

92 


(3) 

(5) 

(’) 


92  COS  ag 

92  cos  Q!5 
<??2  cos  OLi 


(r?j-r,jj)/do 

(^II"  V/'^o 
^^III  ” 


I 


where  ^2^^^  rotation  angle  of  the  i  outer  beam  vibration  at  the  tail  support  and 

(P2  represents  the  vibration  angle  of  the  center  beam  at  branching  point  1.  r}^  is  the 

relative  deflection  of  the  support  point  i  (i  =  I,  H,  III,  IV) .  The  terms  on  the  left  side  of 
equations  (5)  represent  the  relative  outer  beam  rotation  angles  at  the  support  points 
which  affect  the  existing  springs. 


Obviously,  bending  of  the  center  beam  can  cause  axial  forces  in  the  outer  tanks. 
Still,  from  Figure  2,  one  realizes  that  axial  forces  may  attack  the  outer  lox  tanks  only. 
The  outer  fuel  tanks  may  not  resist  axial  forces  because  of  their  special  supports  in  the 
spider  section.  Under  consideration  of  the  axial  lox  tank  forces  S3  Sg  it  follows  from 
Figure  5  that 


Now,  from  Figure  5,  it  may  be  concluded 


(6) 


(92  -  9i)  e  =  S3 
where  e  is  defined  by  Figure  5  and  6^  is  the  lox  tank  influence  coefficient. 
Setting 
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one  obtains 


where 


Equations  3,  4,  8,  9,  and  10  represent  the  complete  boundary  condition  set  of  the 
outer  tanks  at  the  tail  section. 

At  this  point,  some  remarks  will  be  made  about  coupling  between  the  outer  beam 
and  the  center  beam  vibrations.  Obviously,  these  vibrations  can  not  occur  independently 
The  center  beam  distortions  of  branching  point  1  and  2  influence  the  boundary  conditions 
of  the  outer  beams,  while  the  outer  beam  support  forces,  and  moments,  contribute  to 
shear  and  moment  of  the  center  beam  at  branching  points  1  and  2. 
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Coupling  exists  also  between  the  outer  tank  tangential  vibrations.  No  direct 

coupling  exists  between  the  vibrations  of  radial  planes  (1) ,  (2),  (6),  (8).  Hence,  the 

fuel  tank  modes  of  the  planes  (1) ,  (8)  and  the  lox  tank  modes  of  the  planes  (2)  and  f6) 
are  equal  to  each  other.  (See  Figs.  4  and  6) 

In  the  above  boundary  conditions  (i)j  (2),  (3),  (4),  (8)  ,  the  terms  containing 

cos  a.  (i  =  1 ,  2, . . .  3)  and  c^  imply  the  coupling  of  the  outer  beam  vibrations  with  the 
i 

center  beam  vibration.  The  terms  with  the  coefficient  c..,  i,  j  =  3,  5,  7  imply  the 
coupling  of  the  tangential  plane  vibrations  (3) ,  (5) ,  (7),  with  each  other. 

Because 

cos  0^4  =  cos  90®  =  0 

the  vibration  of  the  radial  plane  (4)  is  uncoupled  and  so,  does  not  contribute  to  the  center 
beam  vibration.  Conversely,  a  vibration  in  plane  (4)  can  not  be  excited  since  plane  (4) 
is  perpendicular  to  the  assumed  plane  of  excitation.  For  that  reason,  the  radial  vibration 
of  plane  (4)  could  have  been  omitted  from  this  analysis  without  any  loss  of  generality. 


SECTION  III.  NUMERICAL  RESULTS  AND  COMPARISON  WITH  TEST  RESULTS 


It  is  desirable  to  compare  the  results  of  the  analysis  with  some  dynamic  test 
results  in  order  to  substantiate  the  analysis  procedures  and  the  mathematical  model  to 
the  extent  that  sufficient  structural  detail  was  considered.  Tlie  comparison  can  best  be 
evaluated  if  the  results  of  the  analysis  are  considered  first.  The  analysis,  as  described 
in  this  report,  includes  some  interesting  coupling  of  the  outer  tanks  to  the  main  beam. 

The  dynamic  test  results  which  are  available,  also  need  some  explanation  since  the  test 
was  a  determination  of  transfer  functions  rather  than  a  model  survey. 

Some  typical  mode  shapes  from  the  analysis  are  shown  in  Figures  7  to  16.  (See 
also  Fig.  6)  Each  figure  represents  a  single  mode  of  the  total  vehicle  and  can  be  identi¬ 
fied,  with  some  reservations,  as  either  the  center  line  of  the  vehicle  bending  in  a  first, 
second,  or  higher  beam  mode  shape;  or  the  bending  of  one  of  the  outer  tanks  representing 
the  major  bending.  The  modes  identified  by  vehicle  center  line  are  called  bending 
modes"  and  those  identified  by  outer  tanks  are  called  "cluster  modes.  "  This  identifica¬ 
tion  can  be  rather  useful  in  the  evaluation  of  the  analysis  but  can  also  offer  some  confusion 
in  the  proper  presentation  of  results.  The  identification  is  maintained  because  of  its 
usefulness  in  analysis  and  because  of  the  similar  identification  of  the  dynamic  test 
results. 
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The  complication  of  a  shape  identification  method  can  be  seen  if  the  mode  shapes 
are  studied  as  the  fill  level  is  lowered  in  the  first  stage  of  the  vehicle  to  represent  change 
of  propellant  with  flight  time.  The  simple  and  rather  logical  way  to  show  the  frequency 
trend  throughout  flight  is  to  join  the  lowest  frequency  with  a  smooth  curve,  the  next 
lowest,  and  so  forth  for  the  available  data.  As  the  total  weight  is  reduced,  the  expecta¬ 
tion  would  be  for  an  increase  in  the  frequency  of  the  system.  The  frequency  trend  plotted 
in  Figure  17  shows  that  this  is  true  when  these  next  lowest  frequencies  are  connected. 
When  sufficient  fill  conditions  are  used,  the  frequency  trend  lines  will  represent  the 
frequency  trends  of  the  system. 

Now  the  same  data  points  can  be  joined  by  lines  which  represent  the  same  mode 
shape  trend  from  one  time  point  to  another.  The  trends  of  the  ’’bending’’  modes  and 
’’cluster”  modes  are  shown  in  Figure  18.  The  points  on  these  curves  between  data  points 
do  not  correspond  to  vehicle  frequencies,  with  the  exception  of  the  line  segments  that 
happen  to  be  identical  to  Figure  17.  When  this  identification  is  used,  the  resulting  curves 
can  have  a  decrease  in  frequency  as  the  fill  level  is  lowered.  This  decrease  does  not 
show  in  Figure  18,  but  is  not  unusual.  The  curves  also  intersect  and  cross  each  other  in 
a  manner  which  can  be  described  as  the  cluster  modes  crossing  the  center  line  bending 
modes.  This  crossing  of  the  trends  of  Figure  18  brings  up  the  question  of  crossing  of 
the  trends  of  Figure  17.  These  trends  can  cross,  or  come  together,  since  cases  can 
happen,  and  do  happen  where  two  modes  of  the  system  cannot  be  separated  and  a  single 
frequency  exists  for  two  different  mode  shapes. 

Now  the  test  results  will  be  considered.  The  dynamic  testing  is  conducted  by 
hanging  the  vehicle  on  a  coil  spring-cable  suspension  system  in  a  special  structural  steel 
test  tower.  The  cables  are  attached  to  the  base  of  the  booster  stage  and  to  the  test  tower 
at  about  the  level  of  the  top  of  the  booster.  A  small  angle  is  allowed  between  the  cable 
from  the  vertical,  and  the  springs  between  the  cable  and  tower  are  as  soft  as  practical. 
This  system  provides  little  lateral  restraint  and  only  a  small  end  moment,  which  is 
necessary  to  maintain  the  vertical  stability  of  the  vehicle. 

With  the  vehicle  suspended  in  the  tower,  a  single  shaker  is  used  to  conduct  the 
test  program.  The  shaker  is  attached  at  the  engine  gimbal  station  of  the  first  stage  and 
transfer  functions  are  determined  for  the  one  shaker  location.  The  amplitude  at  the 
frequencies  with  maximum  response  is  plotted  and  decay  damping  at  these  frequencies  is 
determined.  From  these  tests,  the  frequency  trend  curve  and  the  ’’response  mode 
shapes”  will  be  used  for  comparison.  The  comparisons  shown  are  as  complete  as 
practical  and  a  few  comments  are  in  order. 

The  comparison  of  frequency  trends  is  shown  in  Figure  19.  The  first  bending 
mode  frequencies  compare  favorably.  The  largest  deviation  is  in  the  second  mode  and 
can  be  partly  accounted  for  by  the  absence  in  the  analysis  of  an  engine  mode  representation. 
This  has  been  added  to  later  analysis  and  provides  a  better  comparison.  The  test  did 
not  find  an  engine  mode  at  all  test  points  and  is  an  example  of  the  difficulty  in  obtaining 
all  modes  on  the  dynamic  test  vehicle. 


The  frequency  comparison  without  mode  shape  comparison  is  not  sufficient  for 
such  a  system,  therefore,  the  mode  shapes  are  presented  for  the  flight  time  of  35  seconds. 
The  calculated  mode  shapes  show  good  comparison  to  center  line  betiding  of  the  test  modes. 
Some  test  mode  shapes  are  compared  with  two  calculated  modes,  since  the  test  is  a 
"response  mode  shape. "  The  outer  tanks  do  not  compare  as  well  in  amplitude  as  the 
center  line,  since  the  tank  response  is  generally  lower.  The  comparison  to  more  than 
one  calculated  mode  shape  and  the  difference  in  outer  tank  relative  amplitude  can  partly 
be  attributed  to  the  type  of  test  conducted. 

The  results  from  the  analysis  give  good  agreement  with  test  results.  It  can  be 
concluded  that  the  structural  connections  and  general  structural  parameters  have  been 
properly  used  since  the  mode  shapes  compare  favorably  with  test  results.  The  structural 
spring  constants  for  tank  coupling  have  been  improved  since  this  analysis  and  theoretical 
response  analysis  using  the  calculated  mode  shapes  indicate  the  lower  amplitude  response 
of  outer  tanks  is  to  be  expected  and  the  response  of  adjacent  modes  can  be  significant. 


SECTION  rV.  ANALYSIS 


The  transfer  matrix  method  as  outlined  in  the  appendix  will  now  be  applied  to  the 


developed  model.  The  notations  of  the  state  vectors  and  transfer  matrices  are  given  in 
Tables  2  and  3.  The  simple  transfers  of  the  single  beams  (center  beam,  outer  beams) 
may  be  seen  from  Table  2  while  the  transfer  of  the  complete  model  is  given  by  Table  3. 
The  inertia  and  stiffness  matrices  T  and  S  will  be  determined  from  the  givon  vehicle 

.  •  1  I  I  I  (F)  ,  (L) 

data.  (See  appendix.  )  For  this  reason,  the  matrices  Lj,  I-2,  L3,  L2  ,  L2 

which  are  composed  of  inertia  and  stiffness  matrices  (see  appendix)  can  be  considered 

as  known. 


From  Table  3  of  the  transfer  0  —  3  of  the  model  may  be  concluded  as 


where 
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L'l  represents  the  transfer  1^1'  (over  branching  point  1)  and  L'2  represents  the  transfer 
2'  —  2  (over  branching  point  2) . 

To  perform  the  transfer  0  “*3  the  matrices  L'j  and  L2*  must  be  known.  The 
determination  of  these  matrices  represents  the  main  problem  of  the  analysis  and  will  be 
done  in  two  steps.  First,  L'j  will  be  determined  and  then  in  a  direct  way  L2*  .  L'2  will 
not  be  determined  explicitly. 

The  sketch  on  the  top  of  Table  3  shows  the  vehicle  model  represented  by  the 
center  beam  axis.  Because  the  outer  beams  are  removed,  external  forces  and  moments 
must  be  applied  to  ensure  equilibrium. 

It  follows,  then,  from  this  sketch: 


y'l  =  Yi  + 


0 

0 

AMj 

AQi 


(13) 


y2 


y 


0 

0 

AM2 

AQ2 


(14) 


AQ.,  i  =  1,  2  is  the  resultant  of  all  outer  beam  support  forces.  One  may. 


therefore,  write: 


AQ.  =  (-1)^  cos  oil  +  eos 


+  2^  cos  O'.]  ;  i  =  1,  2 

j  =  2  '  J 


where 

;  i  =  1,  2;  j  =  1,  2.  ..  8  are  the  shears  at  the  outer  beam  supports. 
Using  vector  and  matrix  notation,  one  may  write: 


0 

0 

0 

AQ 


-  (-1)  V, 


=  1.  2 


(15) 
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where 


cos  CKi  2  cos  Q'2  .  . 


0  0 

0  0 

0  0 

.  2  cos  a-i  cos  O'; 


;  i  =  1,  2 


AM  is  caused  not  only  by  the  outer  beam  support  reactions,  there  is  also  a  con- 
i 

tribution  of  the  lox  tank  axial  forces.  AM^  can,  therefore,  be  written: 


AM  =  AM!  +  AMV  ;  i  =  1,  2 
ill 


where 


AM'.'  represents  the  moment  contribution  of  the  outer  lox  tank  axial  forces. 

AM'  can  be  expressed  similar  to  equation  15 
i 


(-1)  Vs 


i  =  1,  2 


where 


i  =  1,  2 


and  Vj  results  from  the  matrix  16  by  changing  third  and  fourth  row. 
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The  second  term  of  18  may  be  obtained  from  Figure  5  and  equation  7  as 

1  4-  i  i  +  1  4p^ 

A]Vl”-2(-l)'  "2eS3  =(-!)"  ~  {cp^  -  cp{) 


=  (-1)'^  ^  c  (y2  -  y'i) 


i  =  1,  2 


where 


= 


0  0  0  0 
0  0  0  0 
0  10  0 
0  0  0  0 


In  equation  21  y2  may  be  replaced  by 

f  I  ^ 

Yl  -  1-2  Yi 


(Table  3)  then  it  follows 


=  (-1)^^^  C  A  Yi 


i  =  1,  2 


A  =  u;  (L2  -  U4) 


where  U4  is  the  four  x  four  unit  matrix. 


From  equations  13,  15,  19,  24  it  follows: 


vl  =  Vl  ■  ^4  'll  •  '^3  "’l  +  c  A  Yj 


14 


Now  to  replace  m,  the  boundary  conditions ( 1)  and  (2)  with  regard  to 
j  =  2 . 8,  may  be  introduced.  These  conditions,  written  in  vector  notation  are; 


nij  -  -Kfi+KV2  Yi 


where 


K 


0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 


0  0 
0  0 
K  0 
0  0 
0  0 
0  0 
0  0 
0  0 


0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 


0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
K  0 
0  0 


(26) 


(27) 


f . 
1 


(1)1 


(p. 

1 


h) 


(28) 


1,  2 


V 


2 


0  cos  O'!  0  0 

0  cos  012  0  0 

0  cos  Qig  0  0 

0  cos  0:4  0  0 

0  cos  CKg  0  0 

0  cos  Q!g  0  0 

0  cos  aj  0  0 

0  cos  ttg  0  0 


(29) 


Insertion  of  equation  26  into  equation  25 — and  little  algebra— yields 

(U4  -  c  A)  y;  =  (U4  -  V3  K  V*  )  yi  +  V3  K  fi  -  V4  qj  (30) 


Setting 


(U4  -  c  A) 


=  U 


1  l<2) 

I  -  C  I23 


A  =R 


(31) 
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and 


R(U4- V3  KV2')  =  Bi  1 

RVgKj^-B^  ^ 

R  V4  =  B3  J 

equation  30  changes  over  in 

y  1  =  ^  1  y  1  ■  f  1  ■  ^3  'll 

The  relation  31  may  be  proved  easily  by  verifying 


(32) 


(33) 


(U  -  c  A)  (U  +— ^ - A)  =U 

1  I  ^2) 

1  •  c  '  23 

From  equations  14,  15,  19,  23,  24  one  obtains 

yg  =  (Lj  -  c  A)  yj  +V3m2  +  V4  q2  (34) 


Equation  34  will  be  needed  later  on. 

In  equation  33  the  sixteen  components  of  the  vectors  and  represent  the  un¬ 
knowns  of  the  problem  in  hand.  In  the  following  it  will  be  shown  that  fi  and  can  be 
determined  by  the  outer  beam  boundary  conditions  of  the  tail  section.  From  these  con¬ 
ditions  it  will  be  concluded: 

1  (35) 

'll  =  Q  yi  J 

where  F  and  Q  are  eight  x  four  matrices. 

Then,  combining  equations  33  and  35,  one  obtains 

yj  ^(Bj.B^F-Bg  Q)y4  (36) 

and  hence  (see  Table  3) 

Lj  =Bj  -  B2  F  -  B3  Q  (37) 
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To  determine  the  matrices  F  and  Q  of  equations  35  the  transfers  of  the  outer  beams 
must  be  considered.  In  doing  so,  it  is  convenient  to  use  the  following  notations; 


y,- 


(1) 


y. 


U) 


y. 


(8) 


;  d.2  =fi;  di3  =m^;  di4  - 


i  =  1,  2 


where  q.,  m.,  f.  are  defined  by  equations  17,  20,  and  28,  respectively. 


1  1 


Then  the  outer  beam  transfers  (Table  2)  may  be  written: 


(38) 


d  .  -  y,  L,  -  d  k  =  1,  2,  3,  4 
^2k  ^  ,  k]  1] 

j  =  i 


(39) 


where 


^kj 


kj 

0 

0 

0 

0 

0 

0 

0 


0 


(L) 

kj 


(L) 

kj 


(F) 

kj 


(F) 

kj 


(L) 

kj 


(L) 

kj 


(F) 

kj 


(40) 


k,j  =  1,  2,  3,  4 

Application  of  the  notation  38  to  the  boundary  conditions;  equations  1,  2,  and  26 
respectively,  yields 
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(41) 


di3  -mj  =  -Kf,+KV;'>i 


where  K  and  V2*  are  given  by  equations  27  and  29.  Vi*  follows  from  the  matrix  29  by 
changing  the  first  two  columns. 


Equations  39,  40,  41  result  in 


dji-  =  +  Eifj+L^  qj 

where 

“  Lii  V  ^  +  L 13  K  V2 
^  L  i2  "  Li3  K 


> 


i  =  1,  2,  3,  4J 


(42) 


As  already  mentioned,  the  boundary  conditions  3, 4,8,  will  now  serve  for  deter¬ 
mination  of  fi  and  q^.  Using  the  vector  notation  38  these  conditions  may  be  written: 


•^21  “  V2 

U»(d„.V,%;)=Cj_(y;.y;)tCMd 


23 


where 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

o" 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

;  C  = 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

J-J 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(43) 
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0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

C33 

0 

C35 

0 

0 

0 

0 

1 

0 

0 

0 

0 

C53 

0 

C55 

0 

C57 

0 

0 

0 

0 

1 

0 

0 

0 

0 

Ct5 

0 

C77 

0 

0 

0 

0 

0- 

0 

where 


The  elements  C.  ;  C..;  i,  j  =  3,  5,  7  are  given  by  equations  9  and  iO. 

Lt  Ij 

Insertion  of  equation  23  into  equation  43  and  collection  of  terms  yields 
*^21  “  ^  1  L  2X1 

Us  <^22  “  ^M'^23  “  ^  y  1  \ 

c  =  (U8'v^%Cl)L2-Cl 


Now,  replacement  of  y| ,  d2i,  d22>  ^23  by  equations  33  and  42  changes  the  first  two 
equations  44  over  in 


where 


^21^1^^22^2“^2y2 


^11  ^  14  ^2  ^3 


^12  ^  1  '''  ^1  *-2  ^2 

+vtL2Bi 

N  21  “  Ug'  L  24  +  C  B  3  -  C  ]yiL  34 

N22  "Us  E2'''^^2'  CjyjEg 
R2=-U8D2+CBi+  C]y[D3 
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The  matrices  N..  i,  j  =  1,2  are  eight  x  eight  matrices,  while  R2  are  eight 
X  four  matrices.  Solution  of  this  system  is  given  by  equations  35  where 

F  =  rMRi  -  Nil  N27R2] 

Q  "N27  [-N22  r^Ri  +N22  r^Nii  N27R2  +■'^2]  >  (47) 

I  =  Ni2-Nii  N27N22 
Proof  by  insertion. 

Obviously  the  system  of  equation  45  is  solvable  only  if  the  matrices  I  and  N21 
nonsingular  matrices. 

Now  from  equations  37  and  47  it  follows  that  L;^  is  known. 

To  determine  L2^:'  given  by  the  second  equation  11  and  Table  3,  one  may  start 
from  equation  34.  Considering  equations  33,  35,  and  42  it  follows  from  equation  34: 

y2  =  (J  +  g  F  +  H  Q)  vi  =  L2  vi 
where  J  =  V3D3+V4D4+(L2-cA)Bi 

a  =  ''3E3tV,E,.(q-cA)B^  I 

R  ^34  ^4  ^44  ’  2  ’  ^  ^3 

y 

Because  of  47,  48,  and  49  L2*  is  known  and  hence 
*-  ■  ^3  ^2 

also  (See  eq.  13) .  In  this  way,  the  frequency  equation  of  the  system  can  be  solved  as 
outlined  in  the  appendix.  The  solutions  are  the  eigenvalues 

Xl  =  cuf  ;  A2  =  .  .  . 

Now  the  start  vectors  y  can  be  determined. 

^o 

is  the  start  vector  of  the  center  beam  part  between  top  0  and  branching  point  1. 
The  start  vector  of  the  center  beam  between  branching  points  1  and  2  is  given  by 

LlLiyo 


(48) 


(49) 
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The  start  vectors,  vectors  of  the  outer  beam  vibrations,  are 


The  components  of  these  start  vectors  follow  from  equations  35,  41,  and  47, 

where 

yi  ^  LiYo 
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F - Fuel  Tank 

L - Lox  Tank 


f 


FIGURE  3.  PLANES  OF  EXCITATION 
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Thrust  Outrigger 
Fin  Outrigger 


FIGURE  4.  COORDINATE  SYSTEMS 


Plane  of  Symmetry 


1 


FIGURE  6.  SIGN  CONVENTION  FOR  THE  VEHICLE  MODE  SHAPE  PLOTS 

The  sign  convention  for  the  vehicle  mode  shape  is  positive  as  indicated  by  the 
arrows  in  the  Figure  above.  The  deflection  curves  of  the  tanks  are  plotted  considering 
the  following: 

1.  The  deflection  curves  are  plotted  in  the  planes  as  indicated  above  and  are 
represented  in  the  vehicle  plane  of  symmetry. 

2.  The  reference  axis  for  all  deflections  is  the  undisturbed  main  vehicle  axis. 

3.  The  eight  planes  are  plotted  with  symmetry  considered.  Planes  1,  2,  3,  5, 
and  7  are  plotted  with  their  sign  convention.  Planes  6  and  8  have  the  same  deflections 
as  planes  2  and  1,  respectively,  with  opposite  signs.  Plane  4  has  no  deflection. 

Notice:  This  sign  convention  is  opposite  to  that  used  in  the  analysis. 
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Flight  Time  =  35  sec 


FIGURE  7.  SATURN  I  (SAD-5)  VEHICLE,  RELATIVE  AMPLITUDE  VERSUS  VEHICLE  STATION 


Relative  Amplitude 


FIGURE  8.  SATURN  I  (SAD-5)  VEHICLE,  RELATIVE  AMPLITUDE  VERSUS  VEHICLE  STATION 


Flight  Time  =  35 
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30 


FIGURE  9.  SATURN  I  (SAD-5)  VEHICLE,  RELATIVE  AMPLITUDE  VERSUS  VEHICLE  STATION 


Relative  Amplitude 
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FIGURE  10.  SATURN  I  (SAD-5)  VEHICLE,  RELATIVE  AMPLITUDE  VERSUS  VEHICLE  STATION 


Relative  Amplitude 


Vehicle  Station  (in. ) 

FIGURE  il.  SATURN  I  (SAD-5)  VEHICLE,  RELATIVE  AMPLITUDE  VERSUS  VEHICLE  STATION 


Relative  Amplitude 


FIGURE  12.  SATURN  I  (SAD-5)  VEHICLE,  RELATIVE  AMPLITUDE  VERSUS  VEHICLE  STATION 


Relative  Amplitude 
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FIGURE  14.  SATURN  I  (SAD-5)  VEHICLE,  RELATIVE  AMPLITUDE  VERSUS  VEHICLE  STATION 


Flight  Time  =  35  sec 
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FIGURE  15.  SATURN  I  (SAD-5)  VEHICLE,  RELATIVE  AMPLITUDE  VERSUS  VEHICLE  STATION 


Flight  Time  =  35  sec 


FIGURE  16.  SATURN  I  (SAD-5)  VEHICLE,  RELATIVE  AMPLITUDE  VERSUS  VEHICLE  STATION 


FIGURE  18.  SATURN  I  (SAD-5)  CALCULATED  FREQUENCY  TREND,  MODE  SHAPE 
IDENTIFICATION  VERSUS  FLIGHT  TIME 
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Frequency  Versus  Flight  Time 


—  Test  Frequency 
-  Calculated  Frequency 


NOTES: 


Flight  Time  (sec) 

FIGURE  19.  SATURN  I  (SAD-5)  FREQUENCY  COMPARISON  OF  REPRESENTATIVE 
MODE  SHAPES 
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Matrices 


Table  3.  Transfer  of  the  Center  Beam 
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Table  4.  Used  Matrices  and  Their  Dimensions 
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APPENDIX 


TRANSFER  MATRIX  METHOD 
APPLIED  TO  A  LUMPED  MASS  BEAM 

A  spring-supported  lumped  mass  beam  is  shown  by  Figure  20a.  It  consists  of  masses 

m  and  mass  moments  of  inertia  P.  concentrated  at  Stations  i  on  the  beam  axis  (i  =  1,  2, 
i  1 

3, .  .  .  .  n) .  The  stiffness  between  successive  stations  is  constant. 

If  it  is  assumed  that  this  system  is  free-vibrating,  following  a  sinusoidal  law, 
shear  and  moment  will  obey  this  law  also.  The  equilibrium  conditions  of  the  i^h  mass 
( Fig.  20b)  written  in  the  amplitudes  of  the  mentioned  quantities  are 


M*  =  AP.  q*.  +  M. 
i  111 


Q  *  =-A  m.  y . 
i  11 


(Al) 


where  y^,  9.  are  deflection  and  rotation  angle  at  station  i,  m^,  q..  Q?  represent 

th 

moment  and  shear  at  left  and  right  side,  respectively,  of  the  i  mass  point.  A  is 
the  unknown  parameter,  of  the  problem. 


A  =  0)2  =  4  f2  (A2) 

where  w  is  the  circular  frequency  and  f  is  the  natural  frequency  of  the  vibrating  system. 
The  state  in  a  point  of  the  lumped  mass  system  is  determined  by  the  state  vector. 

y 

9 

^  “  M 

Q 

u  J 
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From  equation  A1  it  follows 


y 

i 

1 

0 

0 

0 

"i 

9.^'^ 

0 

1 

0 

0 

<pi 

1 

0 

AP. 

1 

1 

0 

M. 

1 

Qf!' 

> 

-  Am. 

1 

0 

0 

1 

Qi 

(A3) 


or  y  >:«  =  T  y  (A4) 

i  i  i 

where  the  "inertia"  matrix  is  defined  by  equation  (  A3) . 

A  similar  relation  exists  between  the  state  vectors  y^*  and  y^_^^  (Fig  20c) . 


’'i+i 


(A5) 


where  is  the  stiffness  matrix.  S  ^  can  be  formed  by  solving  the  linear  system  of 
differential  equations: 


GA^.  (y-  -  9)  -  Q  1 

-  El  q> '  =  M  J 


(A6) 


with  the  boundary  conditions 


X  =  X.  ; 
1 


y  _  y.  ;  9=9. 


where  (Fig.  20c) 


M  =  M>i«  +  (x  -  X  )  Q* 
i  i  i 


Q  =  Q.* 

1 


(A7) 


(AS) 
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The  solution  of  equations  A6  (x  =  x.  ^  is  given  by 


/i^  i. 

^i  +  1  "  ^i  ■^i  ‘’’i  "  2EI.  “  (gEI'.  "  GA  . 

1  \  1  SI 

i.  1} 

i.+  l  ^i  El.  i  2EI. 

1  1 


where  i .  =  x.  ,  ,  -  x. 

1  L  +  1  1 


From  equations  A5,  A7,  A8,  and  A9  it  follows 


1  i 
0  1 


0  0 


i?/2EI. 
1  1 

-^./EI. 

1  1 

1 


(i?/6EI. 

i?  /2EI. 

1  1 


i./GA  .) 
1  81 


0  0 


Repeated  application  of  equations  A4  and  A5  results  in 


=  S  T 


1  1  '  1 


=  s.  T.  y.  =  S.  T.  S.  T 


^3  -  "2  '  2  >^2  -  "2  '2 


1  1  '  1 


y.  =  S.  T.  y.  =  S.  T.  ...  S  T  y 

1-1  1-1  'i-i  1-1  1-1  1  1  'i 


^n^n  ^n  ^n-i  ^n-i  ••••  ^i  ^i 


(AlO) 


or  y*  =  L  y 
'n  'i 


where  L  =  T  S  T  ...ST 
n  n-i  n-i  i  l 


(All) 
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The  boundary  conditions  at  Station  1  ( Fig  20a)  are 


Qi  ''  S  yi 


or 


Yl  ~  5i  Qj 


(A12) 


where  k,  k,  and  6,  a  represent  spring  constants  and  influence  coefficients,  respectively, 
of  the  beam  support  at  Station  1. 

Using  matrix  and  vector  notation,  both  types  of  conditions,  from  equation  A12, 
can  be  combined  as  follows; 


t— ^ 

_ i 

1 

o 

o 

1 

I 

1—1 

1 _ 

0  1  -Oi  0 

1 

Ml 

o 

1 

o 

Ml 

Ql 

1 

o 

o 

1 _ 

Ql 

(A13) 


(A14) 


where  the  matrix  E  ,  is  defined  by  equation  A13.  Only  two  of  the  four  columns  of  S  can 
be  applied  —  the  remaining  two  must  be  changed  to  zero  columns.  The  possible  com¬ 
binations  of  applicable  columns  are 


1,  2;  3,  4;  1,  3;  2,  4 


(A15) 


What  combinations  should  be  used  depends  on  the  case  in  question.  Obviously, 
a  free  end  determined  by 


kj  =  Kj  =  0  or  6i  =  oi  -  « 

requires  the  combination  1 ,  2  while  a  fixed  end  given  by 
=  Kj  =  oo  or  6;^  =  cr^  =  0 

the  combination  3,  4  is  required. 
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If  none  of  the  four  quantities  k,  k,  6,  a,  is  zero  or  infinite  any  of  the  above  com¬ 
binations  can  be  selected. 


The  boundary  conditions  of  the  right  end  can  be  written 


+ 


k 

n 


^n 


0 


M 


n 


K 


n 


9 


n 


0 


or 


yn 


+ 


6 

n 


n 


9 


n 


or 

n 


M 


n 


0 


> 


=  0 


Using  vector  and  matrix  notation  equations  A16  change  over  in 


'l 

0 

0 

6  ' 
n 

r - 

o 

\ _ 

0 

1 

-  cr 

n 

0 

e- 

= 

0 

0 

1 

0 

Mn=' 

0 

k 

n 

0 

0 

1 

0 

or 


r 

n  ^n 


0 


(Ale) 


(A17) 


(A18) 


where  is  defined  by  equation  A 17. 

Contrary  to  the  case  of  E^,  two  rows  only  of  E^  can  be  applied  while  the  remain¬ 
ing  rows  must  be  replaced  by  zero  rows.  The  proper  combinations  are  given  again  by 
A15. 


From  the  first  equation  All  and  equations  A14  and  A18  it  follows: 

EnLEjy,=  0  (A19) 

Since  E^  contains  two  zero  rows  while  E^  contains  two  zero  columns,  equation  Al9 
represents  a  linear  homogeneous  system  of  two  equations.  This  system  is  solvable  only 
if  the  determinant  A  equals  zero.  A  is  a  polynominal  in  X.  The  roots  of  A  =  0. 

^2?  ■ 

are  the  eigenvalues  of  the  problem  which  can  be  found  by  a  trial -and-error  procedure. 
Insertion  of  Xi  into  equation  A19  yields  the  solution  y^j  which  represents  the  start 
vector.  Now  from  application  of  equations  AlO  one  obtains  the  state  vectors  of  the 
stations  and  hence  the  i^h  mode  shape. 
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(i)  _ 


^  ^12 


in 


^i2 


9  . 
in 


The  natural  frequency  may  be  calculated  in  accordance  with  equation  A2, 
In  the  case  of  a  free  beam,  equation  A19  may  be  written: 
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The  frequency  equation  then  follows  as: 

^  *^9 


(A20) 


=  0 


(A21) 


Now  X  must  be  determined  so  that  equation  A2l  is  satisfied.  With  this  X,  equation 
A20,  is  solvable  and  yields  the  start  vector. 


‘Pii 

=  0 
lO  - 

The  procedure  to  find  mode  shapes  and  frequencies  is  already  outlined  above. 
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aeronautical  aijd  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.” 

— National  Aeronautics  and  Space  Act  of  1958 
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